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Abstract 

Background: FaWs may occur as unpredictable events or in patterns indicative of potentially modifiable risks and predictive 
of adverse outcomes. Knowing the patterns, risks, and outcomes of falls trajectories may help clinicians plan appropriate 
preventive measures. We hypothesized that clinically distinct trajectories of falls progression, baseline predictors and their 
coincident clinical outcomes could be identified. 

Methods:VJe studied 765 community-dwelling participants in the MOBILIZE Boston Study, who were aged 70 and older and 
followed prospectively for falls over 5 years. Baseline demographic and clinical data were collected by questionnaire and a 
comprehensive clinic examination. Falls, injuries, and hospitalizations were recorded prospectively on daily calendars. 
Group-Based Trajectory Modeling (GBTM) was used to identify trajectories. 

Results: Vile identified 4 distinct trajectories: No Falls (30.1%), Cluster Falls (46.1%), Increasing Falls (5.8%) and Chronic 
Recurring Falls (18.0%). Predictors of Cluster Falls were faster gait speed (OR 1.69 (95CI, 1.50-2.56)) and fall in the past year 
(OR 3.52 (95CI, 2.16-6.34)). Predictors of Increasing Falls were Diabetes Mellitus (OR 4.3 (95CI, 1.4-13.3)) and Cognitive 
Impairment (OR 2.82 (95CI, 1 .34-5.82)). Predictors of Chronic Recurring Falls were multi-morbidity (OR 2.24 (95CI, 1 .60-3.1 6)) 
and fall in the past year (OR 3.82 (95CI, 2.34-6.23)). Symptoms of depression were predictive of all falls trajectories. In the 
Chronic Recurring Falls trajectory group the incidence rate of Hospital visits was 121 (95% CI 63-169) per 1,000 person- 
years; Injurious falls 172 (95% CI 111-237) per 1,000 person-years and Fractures 41 (95% CI 9-78) per 1,000 person-years. 

Conclusions: Falls may occur in clusters over discrete intervals in time, or as chronically increasing or recurring events that 
have a relatively greater risk of adverse outcomes. Patients with multiple falls, multimorbidity, and depressive symptoms 
should be targeted for preventive measures. 
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introduction 

Falls are common among older persons [1,2] and rank among 
the 10 leading causes of death in the United States, resulting in 
more than |19 billion in health care costs annually [3]. Falls 
account for approximately 10% of visits to an emergency 
department and 6% of hospitalizations among Medicare benefi- 
ciaries [4]. 

Although scientific evidence supports associations between a 
number of risk factors and falls [5], efforts to translate these 
findings into effective fall prevention strategies have been limited 
[6]. A Cochrane review showed that multifactorial interventions 



significandy reduce the rate of falls in varying degrees [7,8] but it 
may be difficult to identily the appropriate group to target for 
interventions. It is possible that some people may experience falls 
as random, unpredictable events, while others may have patterns 
of falls that are indicative of potentially modifiable risks factors and 
are predictive of adverse outcomes. Knowing the patterns, risks, 
and outcomes of falls trajectories may help clinicians plan 
appropriate preventive measures. To our knowledge there are 
no longitudinal studies that have observed long-term trajectories of 
falls and determined their predictors and outcomes. 

We hypothesized that: 1- there are distinct clinical patterns of 
falls trajectories ranging from no falls, random fails, clusters of falls. 
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progressively increasing numbers of falls, and chronic recurring 
falls; 2- baseline "predictive factors", such as chronic illnesses 
associated with these patterns could be identified; and 3- clinically 
important health outcomes would be worse in those with 
chronically elevated fall rates. Therefore, we examined a unique 
longitudinal database from the MOBILIZE Boston Study (which 
stands for The Maintenance Of Balance, Independent Living, 
Intellect, and Zest in Elderly (MBS)), which rigorously collected 
falls calendar data from a cohort of community-dwelling elderly 
people over a 5 year period. We used this database to identify 
subgroups of people with distinct falls trajectories, identify baseline 
characteristics associated with these trajectories, and determine 
their coincident clinical outcomes. 

Materials and Methods 

Ethics Statement 

The MOBILIZE Boston Study Trajectory (MBSTraj), Protocol 
Number: 13-029 was reviewed and approved by the Hebrew 
SeniorLife Institutional Review Board (IRB) at the Hebrew 
Rehabilitation Center (HRG) in Boston. Written informed consent 
was obtained from each participant at each phase of MOBILIZE 
Boston Study. The study was conducted according to the 
principles of Helsinki Declaration. 



Study participants 

Study participants were women and men aged 70 years and 
older living in the community in Boston and nearby suburbs. 
Recruitment and enrollment took place from September 2005 to 
January 2008 within a defined geographic area bounded by a 5- 
mile radius around the Institute for Aging Research at HRC in 
Boston. The sampling area was chosen to capture a diverse urban 
and suburban population, to increase likelihood of recognition of 
the study center, and to minimize transportation burden. Details 
of the study methods were published previously [9,10]. Initial 
eligibility was based on age 70 years or older, ability to walk 20 feet 
without personal assistance, ability to communicate in English, 
and the expectation of staying in the area for 2 years. Following 
the initial recruitment visit, study staff contacted prospective 
enroUees by telephone to confirm eligibility and schedule the 
baseline home and clinic visits. During the home visit, written 
informed consent was obtained and participants were screened 
and excluded for moderate or severe cognitive impairment using 
the Mini-Mental State Examination (MMSE score, <18) [11,12]. 

Falls and Clinical Outcomes Assessments 

During the Follow-up, a fall was defined as unintentionally 
coming to rest on the ground or other lower level not as a result of 
a major intrinsic event (e.g., myocardial infarction, stroke, or 
seizure) or an overwhelming external hazard (e.g., hit by a vehicle) 
[13]. Participants were instructed to complete and return monthly 
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Figure 1. Patterns of Falls Trajectories over 5 years in Older Adults (The MOBILIZE Boston Study). — Solid lines are estimated 
trajectories. - — Dashed lines are observed trajectories for each group. 
doi:1 0.1 371/journal.pone.01 06363.g001 
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Table 1. Summary of Baseline Characteristics (The MOBILIZE Boston Study). 





Demographic Characteristics and Risk Factors 


Total sar 


nple n = 765° 




No. 


% 


Age, mean (SD), y 


78.1 


(5.4) 


Women, 


488 


63.8 


Race, 


White 


609 


79.6 


Nonwhite 


156 


20.4 


Years of Education, mean (SD), y 


14.8 


(6.1) 


Body mass index, l<g/m^, ^ 


<25 


233 


30.4 


25-29.9 


297 


38.8 


2:30 


206 


26.9 


Fell in past year. 


289 


37.8 


Comorbidities 2:2, 


491 


64.2 


Hypertension, 


607 


79.3 


Previous stroke, 


72 


9.4 


Diabetes mellitus. 


172 


22.5 


Hyperlipidemia, 


469 


61.3 


Congestive heart failure, 


38 


5.0 


Mini-Mental State Examination score <24, 


154 


20.1 


CES-D Revised score, mean (SD), points 


11.0 


(11.2) 


Physical activity score, 


0-66 


245 


32.0 


66.01-124 


262 


34.3 


124.01-559 


258 


33.7 


Impaired balance, score <4 out 7, ' 


318 


41.6 


Slow gait speed, <0.78 m/s, ^ 


189 


24.7 


Slow chair stands, > 16.37 s, ^ 


95 


12.4 



^Forty four of the 809 older adults were missing data and therefore excluded from the analysis. 
^Body mass index is calculated as weight in kilograms divided by height in meters squared. 
^Mini-Mental State Examination (MMSE) cut off point for cognitive impairment [11,12]. 
^Physical activity tertiles measured using the Physical Activity Scale for the Elderly [15]. 

^Balance score was based on 4 progressively difficult stands: feet side by side, semi-tandem, tandem, and 1-leg stand [17]. 
^Slow gait speed (m/s) is the lowest 25% based on time of lowest of 2 usual-paced 4-meter walks [17]. 
^Highest quartile (slowest performance) of time to complete 5 repeated chair stands [17]. 
doi:10.1371/journal.pone.0106363.t001 



falls calendar postcards designed to be posted on a refrigerator. On 
the postcards, participants were to record an F for each fall on the 
day it occurred and an A'^ on days when No Falls occurred. This 
approach has been well-validated for use in epidemiological cohort 
studies [14]. The calendar postcards also included questions about 
whether the participant experienced a hospitalization during the 
preceding month. Research staff monitored the return of the 
calendars and on any given month, approximately one-third of 
participants were called for missing or incomplete calendars. AU 
subjects who reported falls were also called to determine the 
circumstances of the fall and clinical outcomes whether any 
injuries (e.g. fractures) and hospital visits were incurred. 

Covariates 

Covariates included sociodemographic characteristics, physio- 
logic risk factors, health status, and amount of physical activity. 
Sociodemographic characteristics assessed in the home interview 
included age, sex, race (self-identified), and years of education. 



Cognitive status was assessed using the MMSE, scored 0-30 [12]. 
We used the validated Physical Activity Scale for the Elderly 
(PASE) to measure physical activity in the previous week [15]. 
Participants were asked about physician-diagnosed major medical 
conditions. Details of the study methods were published previously 
[9,10]. During the clinic visit, several measures were calculated. 
Diabetes was defined using an algorithm based on self-reported 
diabetes, use of antidiabetic medications, and laboratory measures 
from the baseline clinic visit including random glucose (200 mg/ 
dL) and hemoglobin Ajj. (7%). Depression was assessed using the 
Eaton method based on a modification of the 20-item Centers for 
Epidemiologic Studies Depression scale [16]. Body mass index 
(calculated as weight in kilograms divided by height in meters 
squared) was calculated from measured height and weight. 
Standing balance was scored using 4 timed tests (side-by-side, 
semi-tandem, tandem, and 1-leg stands) [17]. For the timed chair 
stands test, participants were asked to fold their arms across their 
chest and stand up and down from a chair 5 times as quickly as 
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possible [17]. Gait speed was based on the shortest time of 2 trials 
of a usual-paced 4-meter walk [17]. 

Data Analysis 

To identify chnically distinct trajectories of falls, we use Group- 
Based Trajectory Modeling (GBTM) [18] . This method allowed us 
to simultaneously estimate probabilities for multiple trajectories 
rather than a single mean for the population, as is the case for 
traditional regression or growth-curve models. We used SAS 
software, and the PROC TRAJ macro (http:/ /www.andrew.cmu. 
edu/user/bjones) [19,20], a closed-source module developed 
specifically for use with SAS software, which fits a semiparametric 
mixture model to longitudinal data with the use of the maximum- 
likelihood method. It is possible that some people may experience 
falls as random, unpredictable events, while others may have 
patterns of falls that are indicative of potentially modifiable risks 
factors and are predictive of adverse outcomes. Knowing the 
patterns, risks, and outcomes of falls trajectories may help 
clinicians plan appropriate preventive measures. The metric for 
defining the trajectory was years into the study. The number of 
falls that each participant reported each month was summated for 
each year over the 5 years study period and quantified as a 'Tails 
Per Year of Aging" Score (FPYA Score). For each one year time 
interval. No Falls equaled FPYA 0, one fall equaled FPYA 1, two 
falls equaled FPYA 2 (recurrent falls), and three or more falls 
equaled FPYA 3 (high number of recurrent falls). FPYA scores 
were examined as a function of different time periods from 1 to 5 
years long, independent of the chronological time at which these 
periods started (e.g., a 3 year period could start at year 1 and 
include years 1-3, year 2 and include years 2-4, or year 3 and 
include years 3-5 of the study). PROC TRAJ was used with the 
follow-up time metric from 1 to 5 years long independent of their 
age. So interval time was 1-year and Falls Per Year Aging (FPYA) 
was scored with 1-year follow-up time metric. FPYA (Falls per year 
aging) as the y-axis in Figure 1 represents the average of fall rate 
per year as a function of different intervals in time (shown on the x- 
axis). These scores were modeled as a censored normal 
distribution. We used the Bayesian Information Criterion (BIC) 
to test from three to six trajectories and to determine whether each 
trajectory was best fit by intercept only (i.e., constant) or by linear, 
quadratic, or cubic terms [2 1] . The final model was selected based 
on a combination of the Bayesian information criterion (BIC; 
where the value closest to 0 indicates the best-fitting model) and by 
estimated trajectory group proportions that were sufficiently large 
(e.g., 0.05) [22]. These analyses were repeated after adjustment for 
age, sex, race or ethnic group, years of education, chronic 
conditions such as hypertension, stroke, diabetes mellitus, hyper- 
lipidemia, cognitive status and previous falls at baseline using the 
Proc Traj software "risk" command. The proportions of older 
adults classified according to each trajectory, the mean probability 
of assignment, and the proportions with poor fit were based on the 
original data, and 95% confidence intervals were estimated with 
the use of 1000 bootstrap samples [18]. Missing data on the 
trajectory modeling and on time-varying variables were handled 
with a maximum likelihood approach [18] together with the 
missing-at-random assumption, which assumes that for each 
individual, the likely values for missing data on the trajectory 
and time-varying variables can be estimated from other available 
observed data. This approach uses available information on each 
case for constructing the trajectories rather than deleting 
individuals with missing observations. 

After determining the most appropriate GBTM, group status 
for each individual was obtained to identify relevant predictors of 
each falls trajectory group. An individual was assigned to the 
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Figure 2. Clinical Outcomes according to Falls Trajectory group membersKiip (TKie MOBILIZE Boston Study). A. Incidence Rate of 
Clinical Outcomes according to Falls Trajectory Group. B. Box plot showing Averages of Clinical Outcomes according to Falls Trajectory Group: 
Trajectories: 1 = No Falls; 2 = Cluster Falls; 3 = Increasing Falls; 4= Chronic Recurring Falls. Injurious Falls: 2&4 p<0.001 (S); 3&4 p<0.001 (S); All others 
NS. Hospital visits: All are NS. Fractures: 1&2 p<0.001 (S), 1&3 p<0.001 (S), 1&4<0.0001, All others are NS. 
doi:1 0.1 371/journal.pone.01 06363.g002 



trajectory group in which he or she was most likely to be, as 
determined by the group posterior probabilities from the final 
model. After using the "risk" command in Proc Traj software to 
identify predictors with beta coefficients and p values, we 
confirmed the findings and computed odds ratios using multino- 
mial analyses. The outcome variable was falls trajectory group 
membership, and the testable predictors were sociodemographic 
characteristics, physiologic risk factors, health status, and amount 
of physical activity. The models were estimated using PROC 
LOGISTIC in SAS version 9.3. Partial was calculated for 
statistically significant predictors using R^max values [23,24] and 
the Area Under the Curve was calculated for the model 
discrimination [25]. We assessed relevant clinical outcomes such 
as injurious falls, falls resulting in fractures, and hospital visits 
according to falls trajectories. All analyses were performed using 
SAS software, version 9.3 (SAS Institute Inc, Cary, North 
Carolina). A two-sided P value of less than 0.05 was considered 
indicative of statistical significance. 

Results 

Participants 

Eight hundred nine individuals met study inclusion. Forty four 
of them older adults were missing data and therefore excluded 
from the analysis. Table 1 summarizes the demographic and 
health characteristics of the study sample at baseline. Of the 765 
older adults included in this study, 276 (36. 1 %) were male and 489 
(63.9%) female, 156 (20.4%) were nonwhite, 245 (32.0%) met the 
PASE score of low intensity physical activity. On average, the 
cohort was 78.1 ±5.4 of age (range 64—97). Over five years, 90 
(11.8%) older adults died during the follow-up period, with an 
incidence rate of 24 (95% CI 13-34) per 1,000 person-years. 

Determining Falls Trajectories over 5 years 

The BIC values and estimated proportions for the 3-, 4-, and 5- 
group GBTM models shown in Table 2 were used to determine 
the best fit. A 6-group model was tested but failed to converge. 



Consequentiy, the 4 group model was selected for each fall 
trajectory group. Additional diagnostic criteria for judging the 
adequacy of a GBTM demonstrated that the 4 group model 
performed well based on the Nagin criteria [18]. Figure 1 
illustrates each of the 4 trajectories along with the average raw 
group data at each point. The 4 distinct trajectories were: Cluster 
FaUs (46.1%; 95CI, 40.0%-48.8%; n= 353) P<0.001, Increasing 
FaUs (5.8%; 95CI, 3.7%-11.4%; n = 44) P<0.001, Chronic 
Recurring Falls (18.0%; 95CI, 12.7%-19.1%; n= 138) P<0.001, 
and No Falls (30.1%; 95CI, 28.8%-35.6%; n = 230) P<0.001. 
Table 3 presents descriptive data for each falls trajectory group. 

Coincident Clinical Outcomes of Falls Trajectories 

Figure 2 panel A shows the distribution of the incidence of 
clinical outcomes namely injurious falls, fractures and hospital 
visits for each falls trajectory and summarized below. The average 
of clinical outcomes and incidence followed the same distribution 
in falls trajectories. The incidence oi Injurious FaUs was 172 (95% 
CI 111-237) per 1,000 person-years in Chronic Recurring Falls 
trajectory group, 126 (95% CI 98-169) per 1,000 person-years in 
Cluster Falls trajectory group, and 1 12 (95% CI 51-267) per 1,000 
person-years in Increasing Falls trajectory group. The incidence of 
Fractures was 41 (95% CI 9-78) per 1,000 person-years in the 
Chronic Recurring FaUs trajectory group, 21 (95% CI 7-38) per 
1,000 person-years in the Cluster FaUs trajectory group, 23 (95% 
CI 10-67) per 1,000 person-years in Increasing FaUs trajectory 
group, and 2 per 1,000 person-years in No FaUs trajectory group. 
The incidence oi Hospital visits was 121 (95% CI 63-169) per 
1,000 person-years in Chronic Recurring FaUs trajectory group, 
116 (95% CI 20-207) per 1,000 person-years in Increasing FaUs 
trajectory group, 99 (95% CI 68-130) per 1,000 person-years in 
Cluster Falls trajectory group and 97 (95% CI 58-134) per 1,000 
person-years in the No FaUs trajectory group. Figure 2 panel B 
shows that during the 5-year follow-up, the mean number of 
injurious faUs was significantiy higher in the "Chronic Recurring 
FaUs" trajectory group. The mean number of fractures was also 
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higher in the "Increasing Falls" and "Chronic Recurring Falls" 
trajectory groups compare to the "No falls Group". 

Predictors of Falls Trajectories 

Multinomial logistic regression revealed that specific baseline 
characteristics predicted membership within each of the three 
"faUer" trajectory groups as compared to the no falls trajectory 
group are shown in Table 4. Predictors of the Cluster Falls 

trajectory were faster gait speed and falls in past year. Predictors of 
the Increasing Falls trajectory were Diabetes MeUitus (DM) and 
Cognitive Impairment. Predictors of Chronic Recurring Falls were 
two or more comorbid chronic conditions and falls in the past 
year. Symptoms of depression were predictive of all groups. 

Discussion 

In this 5-year prospective cohort study, we examined the course 
of falls in community-dwelling older adults. We found four 
clinically distinct trajectories of falls that we have labeled: cluster 
falls, increasing falls, chronic recurring falls, and no falls. We found 
that Diabetes mellitus, cognitive impairment, fast gait speed, falls 
in the past year at baseline, and multi-morbidity are predictors of 
these trajectories. Importantly, these trajectories were associated 
with adverse outcomes including injurious falls and fractures. 
Notably, people with chronic recurrent falls had the highest rate of 
injurious falls and fractures. There was a marginally significant 
increase in hospitalizations among those with increasing falls. 

The number of falls that each participant reported each month 
was summated for each year over the 5 years study period and 
quantified as a "Falls Per Year of Aging" Score (FPYA Score). For 
each one year time interval. No Falls equaled FPYA 0, one fall 
equaled FPYA 1, two falls equaled FPYA 2 (recurrent falls), and 
three or more falls equaled FPYA 3 (high number of recurrent 
falls). These scores were distributed according the censored normal 
distribution and their modeling identified 4-groups. 

The major utility of the 4-group solution is in the ability of these 
4 groups to predict adverse outcomes. In addition, the 4-group 
trajectory solution identified some risk factors to target for 
interventions. These patterns of falls had a few distinguishing risk 
factors described above. Some of these are potentially modifiable, 
such as diabetes for increasing falls and comorbid conditions for 
Chronic recurring falls. The greater proportion of fast walkers in 
the cluster group may represent people who are less careful or 
have more acute condition rather than chronic diseases (e.g. these 
comorbit)' scores are lower than other fall groups). We believe that 
falls incidence declines in the chronic recurrent falls group 
beginning in year 4 because people who had recurrent falls had 
chnically adverse events like hospital visits or fractures which 
reduced their mobiUty and chance of falling. Recurrent falls over 
3-years may signal the begirming of fLmctional decline. 

This study is unique in demonstrating that not all falls have the 
same clinical implications. As we have shown, some may occur in 
relatively short-lived clusters, perhaps from an acute illness or 
exposure to a new drug or environmental hazard; some may 
increase in frequency over time, possibly due to the progression of 
a disease or disability; and some may recur at a relatively high rate 
over time. Each of these presentations appears to have a distinct 
set of risk factors and in some circumstances different prognostic 
implications. For example, chronic recurrent falls seem to occur in 
people with multiple co-morbidities (multi-morbidity) [26,27] and 



depressive symptoms, and, as might be expected, are associated 
with a relatively high rate of injuries and fractures. 

Although a high incidence of falls has been reported in older 
populations [1], little is known about the progression of falls over a 
long period of time. To our knowledge there are no longitudinal 
studies that observed the trajectories of fallers. While some falls are 
isolated, non-recurring events, others may mark the beginning of a 
progressive down-hill course in which the injuries and fear 
associated with an initial fall may precipitate recurrences and 
spiraling functional decline. A previous study showed that illnesses 
and injuries leading to either hospitalization or restricted activity 
are strongly associated with the initial onset of disability [28]. 
Another study extended this earlier work by demonstrating that 
exposure to intervening illnesses and injuries are also associated 
with the subsequent course of disability [29]. Therefore, it is 
important to distinguish those falls that are likely to be isolated 
episodes from those that lead to further disability. A Cochrane 
review highlighted the difficulty identifying an appropriate group 
for preventive efforts as a result of this heterogeneity in the 
presentation and subsequent course of falls in older populations 
[7,8]. 

Our study has several limitations. 1) We used baseline clinical 
characteristics that were obtained up to .5 years before the falls 
occurred, so we were unable to identify risk factors at the time the 
falls. This limited our ability to determine whether cluster falls 
were due to acute illnesses, drugs, or environmental hazards, and 
whether other patterns followed the course of chronic diseases. 2) 
Some of our participants experienced falls before the study began. 
Thus, these falls were left censored in our analysis. It is therefore 
possible that some of our non-fallers experienced fall episodes 
before the study began. We performed secondary analyses, 
stratifying previous fallers by non-fallers and found similar 
patterns. Unfortunately, the sample sizes of these subgroups were 
too small to draw meaningful conclusions. 3) Finally, because our 
study participants were only from a Massachusetts geographic area 
within a 5-mile radius of the Institute for Aging Research, the 
cohort may not be generalizablc. However, a comparison of the 
demographics of MBS study participants to the US Census showed 
comparable distributions by sex and racial group in the population 
aged 70 and older [9]. 

In summary, our study provides evidence suggesting four 
clini[:ally distinct patterns of falls over fi\c 5'cars among 
community-dwelling older adults. Each of these patterns has a 
distinct set of independent risk factors. Compared to all other 
groups, chronic recurring falls are associated with the highest rate 
of injuries. Older adults with recurrent falls should be targeted by 
health care professionals to identify their underlying cause and 
implement interventions to prevent subsequent injury. 
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